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Abstract 


This essay considers an optimal control problem (OCP) governed by a 
system of Fredholm integral equations (FIE). In this paper, collocation 
approach with utilizing Lagrange polynomials is introduced to transform 
the OCP into a nonlinear programming problem (NLP). An efficient op- 
timization method in Mathematica software is utilized to solve NLP. The 
convergence analysis is discussed, which show the theoretical structure 
behind the propounded technique under some assumptions. In this es- 
say, computational outcomes are given to demonstrate the adaptability, 
forthrightness, and relationship of the calculations manufactured. A prac- 
tical real-world problem involving hanging chain in classical mechanic is 
also dissolved utilizing the approach proposed. 
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1 Introduction 


The mathematical formulations of most physical and mechanical engineering 
models lead to FIE. FIEs occur in different real world models like the mass 
dispensation of polymers in polymeric melt, linear forward modeling, signal 
processes etc [7, 26]. The OCPs inclusive integral equations are a different 
class of OCPs that have a substantial role in population dynamics, econ- 
omy, heat conduction, continuum mechanics, electrical power centers, water 
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resource management, and mass transport [13, 19]. In particular, optimal 
control of systems managed by FIE is significant in usages like the OCP 
regarding the Ornstein-Uhlenbeck process which arises from statistical com- 
munication theory [18]. In recent years, many methods for solving OCPs of 
FIE, whether direct or indirect, have been proposed. One can read about 
this subject in references [2, 1, 8, 10, 16, 20, 24]. 

ProblemB: Determine the optimal control vector u* and the correspond- 
ing optimal state vector x* for a category of OCPs governed by a system of 
nonlinear FIEs described as follows: 


T 
min J(z,w) =| F(t,x1(),-++ ,@n(t), ur(t),- ++ , Un (t)) dt, (1) 


subject to 


T 
x;(t) = y(t | K;(t, 8,21(s),-++ ,@n(s), ur(s),--+ ,Um(s))ds, jg =1,--- ,n, 
(2) 


where x = [21,22,::-,2n]7 € RK", wu = [ui,ue,::: um]? € R™, F € 
C((0, T], R", RK”) and K; € C((0, T], [0, 7], KR”, RK”). 

Many researchers have been solved and carefully examined non-linear system 
of FIEs by utilizing orthogonal basis functions. For example, see the refer- 
ences [3, 4, 21, 22]. The orthogonal nature causes the solution of the prelimi- 
nary problem to transmute into an algebraic system of equations. More ever, 
by utilizing an operational matrix of integral or derivative of these functions, 
the computational convolution of the problem is substantially diminished. 


For dissolving the mentioned problem, the Lagrange polynomials and 
Gauss—-Legendre (GL) quadrature rule are used. A robust direct approach 
is established to transform the problem 6 into a NLP. This method is based 
on the expansion of control and state functions with Lagrange polynomials 
together with the GL integration method. We have utilized Wolfram Math- 
ematica 12 for obtaining the solution of the NLP. Convergence analysis and 
related theorems are discussed in details. The reliability of the approach is 
investigated numerically by dissolving some illustrative examples. The OCP, 
which modeled the problem of classical mechanics for power lines are also 
described and solved by propounded approach. 


The remaining structure of our article is given as follows: In Section 
2, by using approximation of control and state variables, we obtain a NLP 
corresponding to problem B and also discuss about the method of solving 
the resulted NLP. Convergence analysis is given in Section 3. In Section 4, 
we demonstrate the efficiency of the propounded method to several examples 
and show the computational results of the method. Finally, in Section 5, 
concluding comments and future extensions are summarized. 
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2 Proposed method 


Consider the problem 6 in equations (1) and (2). In this section, by utilizing 
the properties of Lagrange polynomials and GL integration, the considered 
OCP is changed to a NLP, then the function NMinimze in Wolfram Mathe- 
matica software is used to solve the resulted NLP. 


2.1 Discretization of problem B 


Firstly, we approximate the control function u(t) and corresponding state 
function x(t) over the interval [0,7] by the following linear combination of 
the Lagrange polynomials 


N 
x(t) od £5 y(t) od S > 2;,Qi(t) = X;" Q(t), J=Hl,---,n, (3) 
1=0 


and 


N 
ux(t) ~ Urn (t) = S- uniQi(t) = Ux’ Q(t), k=1,---,m, (4) 
1=0 


where Q(t) = [Qo(t), Qi(t),--- ,Qn(@)|?, X} = [xj0, 251, 32,+°+ , ey]? and 
Ux = [uxo,Uk1,UK25°°* »UkN]?, in which xj, and ux; for i = 0,1,---,N, 
jg =1,--:,n and k = 1,--- ,m are the unknown coefficients and, Q;(t) are 
the Lagrange polynomials which are defined as follows: 


7 t-t, w(t) 
aw= I] t;-t, (t—t)w'(ti)’ 
O0<s<N 


where w(t) = ll ve —t,;) and t;, for i = 0,1,2,---,N, are the zeros of 
Legendre polynomials [11] transformed to the interval [0,7]. The transfer is 
done according to the interval in which the problem is defined. By discritizing 
the FIE given in Eq (2) by using (3) and (4) and collocation nodes t, ,we 
obtain for 7 =1,--- ,nandq=0,---,N 


a 
rjn (tq) = y(tq) +f Kj (tq, 5, r1n(s), a »Tnn(8), uin(s), oa ,Umn(s))ds. 
(6) 


For using GL quadrature formula, a simple linear transmutation must be 
produced with the following form 


ee Sut), (7) 
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Then, the dynamical system (2) is converted to 


7 T y 7 
£5N (tq) = y(tq) ai 9 / i Kj (te, Pate ae 1), win(>(H + 1), ie tan (> (u + 1) 
T T 
un (SUF Da tha (4+ 1) de 
(8) 
We apply the GL quadrature for approximating the integral given in Eq (8) 
ie 7 i 7 
3 (tq) = y(tq) a 9 So wk; (tq, 5 (ma ole 1), tin (5 (Ma + 1), - nn (> (MI + 1)); 
1=0 
T ah 
win (> (Ht + 1)),-++ Um (G (Ma + 1))), 
(9) 


where wy is defined as follows: 


2 
(ie (42)?) (Pag (un) 
in which Py +1 is the Legendre polynomial of order N +1 [11]. By using GL 


quadrature and equations (3) and (4), we can approximate the cost functional 
in equation (1) as follows: 


W= 


joa [reno era st (f)) db = 5 
[Ff (Gus 1), 21 (Gus 1) st? 5 Ca (Gus 1) ,ULN (Gus 1) 
oti (U1) de oe (Flot Dew (For +2), 


tan (Fa + 0) stay (S42) s+ stim (Fe +29) ). 


Eventually, the OCP governed by a system of FIE in problem 6 has been 
reduced to a NLP by (10) as objective functional and (9) as constraints. We 
call this NLP as problem By. 
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2.2 Optimization algorithm 


In this paper, for solving NLP, which we call as problem By, given in equa- 
tions (9) and (10), the Nminimize function in Wolfram Mathematica has 
been utilized. Nminimize find global optimum numerically. This function 
implement several algorithms for solving optimization problems according to 
the kind of optimization problem. Now by solving this NLP and determining 
the vectors X; and U; for 7 = 1,---,n and k = 1,--- ,m, we can find a 
numerical solution of problem B given in equations (1) and (2). 


3 Convergence analysis 


We suppose that the finite set, Py, consists of all Lagrange polynomials from 
degree less than or equal to N. The set of Py is linearly independent and 
a basis of the vector space for polynomials of degree at most equal to N. If 
F(t) is a function in L?[0,T], F(t) has a best estimation Fo € Pn, so we 
possess [17] 

VgeELl?0,T): ||F-Folle < IlF — alle. (11) 


Suppose S,, € Py, then there is coefficients s;, 1 = 0,1,--- ,N, whereas 
N 
Sn = )° 81 Qi(t) = S* Q(t), (12) 
i=0 


where S = [59,51,°:- , Sn], and Q(t) = [Qo(t), Qi(t),--- , Qn(d)]. 


Theorem 1. Let F € L?(0, 1] be estimated by Lagrange polynomials {Q;}/V, 
that is, Fy := > C;Q;(t). If ex is absolutely error of Fy, then ||F — 
Fy\lo < (Way and limy_,.o en = 0 and by increasing N; the errors quickly 
tends to zero. D is a constant. 


Proof. The proof is given in [23]. 


Theorem 1 indicates the accuracy estimation of Lagrange polynomials. 


Definition 1. The vectors x(t) € R” and u(t) € R”™ are called admissible, 
if they satisfy in equation (2). The set of admissible pairs {(x,u)}, for j = 
1,---,n, is defined as follows: 


Tv 
wy = {(x(t), u(t))| 25 (t) = yy) +f Kj(t,8,21(8),-+-, 
Ln(S),U1(S),-++ ,Um(s))ds}. (13) 
Define U as the admissible control functions. Let 


EN = (Lin, 22N,°*',%nn), UN = (UN, U2N,*** ,UmN); (14) 
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then we have 


and 


T 
dy = {(en(t),un(®)|eiw(t) = u(t) +f Kilt,s.2n(9) 
+ @nv(s), tin (8),+** stima(s))ds}. (16) 


The proof of the following theorem is given by utilizing of the same argument 
in [25], but some revisions are given according to the problem of our paper. 


Theorem 2. Assume J; = infy, J and J* = infy J and J* is finite and 
unique, then the following relation is hold. 


dy de = ed) Sa So =i ad), (17) 


Proof. From the definition of ~;, the following relation is hold 


Wy Coe C307? Cn Ce Cy. (18) 


According to (18), {J} is a non-increasing sequence that is bounded from 
below, so it is convergent to J*. 


Remark 1. The computational distance can be transmuted from [0,T] to 
any ideal interval via an affine alteration. 


Lemma 1. Suppose that ¢; and w; for 1 <7 < N are GL nodes and weights 
respectively, Let f(t) be Rimann integrable. Then [9] 


N 


[soa = lim S° f(ti)wi. (19) 


N-co ¢ 
i=1 


Theorem 3. Assume that 7* € #” is exact state solution of equation (2) 
with u* € #™” as exact control functions, let vi, and uy, be numerical solu- 
tions of 


T 
LN (tq) = y(tq) +f K (tg, 8,®1n(S) +++ ,&nn(S),Uin,++: ,Umn))ds. (20) 


Also assume that the function K(t, s,21(s),--+ ,&n(8),U1(S),°++ , Um(s)) sat- 
isfies Lipschitz condition with respect to variables x and wu as follows: 


|K(é, s, x1(s), a ;in(s), Ur(s), a , Um(8)) 


~ K(t, s,2'1(8), aS ,2'n(s),u'1(s), bese ,U'm(8))||2 
< Lx{\|21(s) — 2'1(s)|l2 +--+ + |lan(s) — 2’n(s)|l2 + [lua (s) — w'1(s)Il2 
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+++ + ||Un(s) — u'm(s)]l2}, (21) 


where Lx is Lipschitz constant. Then we have: 


le) —exltnlla s (Ga) (22) 


where A is a constant which is non-dependent to N. 


Proof. We have for 7 = 1,:--,n 
T 
let(t) — aty()II < if K (t,s,27(s),--- ,2%(t),ut(s),--- ,u%,(s)) ds 
T 
=a K (t,8,2ty(s),--+ -2%y(8),Uty(8),-++ -whv(s)) dll2 
T 
< f[K (t.s,2%(s)oo+ .25(0), wi(s)e--> uh (9)) ds 


0 
=e (t, 8, rin(s), a Cnn (8), Uin(S), a ,Umn(8)) ds||2 
T 
< ix f ||zi(s) — tin (s)|l2 +--+ + |len(s) — rn (s)Il2 
+ ||ui(s) — uin(s)|l2 +--+, [lum (s) — umn (s)|leds 


nm 97 
=> | \l22(s) — aty(s)]lods 
i=l 
m aT 
+ ie | ful (9) — ubyy(2) lade 
i=1 79 


(23) 


< LxT) | |le?(s) — chy(s)|l2 + LxT >| |lui(s) — uiv(s)lle 


(24) 
ei Lj A L; 
< 2LKT (aa) +2LKT (524) 
<n (S34 +n) < arene 
(25) 


In (25) A = 2LeT (Dy Lit+ DM Li), so by increasing N for j = 
1,2,-++ ,n, ||x}(t) — 2} y(t)|| goes to zero. 
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4 Numerical experiments 


In this section, the proficiency of the proposed technique is shown by applying 
to some given numerical examples of OCP for systems governed by FIE. 


Example 1. Observe the following OCP 
1 
ion TG aN = i (x(t) — 0.8182 — 2.727342)? + (u(t) — #)2dt, (26) 
0 


subject to 
1 
z(t) = t? +f (t? + u(s))a(s)ds, (27) 
0 
with the exact solutions, u*(t) = t? and x*(t) = 0.8182 + 2.7273¢?. 


By exerting the propounding method, the approximate solutions are 
approximately x*(t) = 0.8182 — 8.78094 « 10~!°t + 2.7273t? and u*(t) = 
1.47439 * 10-6 + 1.7. The optimal value of J* for N = 2 is 1.20133 * 1071. 


Example 2. Investigate the following minimization problem 


J (#1 (t), £2(t), u(t)) = i; (ay (t)—t? +1)? + (xo(t)—14t?)?+(u(t)—sin(t))7dE, 


-1 
(28) 
subject to FIE system 


a(t) = —1+4+ / t?u(s) (a1(s) + 22(s)) ds, 


1 
#9(t) =1-24 ; (s — u(s))(a1(s) + 22(s))ds. 
-1 
The precise optimal solutions of this OCP are x(t) = t? — 1, 23(t) =1-¢? 
and u*(t) = sin(t). 

After dissolving with the propounded approach, we obtain the following 
result 


gy*(t) = —1+5.55112 x 10-¢+ 2, 


wo*(t) = 1+5.55112 x 1071" — 12? 
u*(t) = —1.50469 x 10~*4 + 0.902957¢ — 1.11022 x 1071%¢?, 


which approximates the exact functions with good approximation. The 
Maclaurin series of sine function is sin(x) = # — ae + =r —**>« AS We can see 
in the approximate control function u*(t) = —1.50469 x 10724 + 0.902957t — 
1.11022 x 107!°t?, the constant sentence and the coefficient of sentence in- 


cludes x? is approximately equal to zero. The optimal value of objective 
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function for different value of N is given in Table 1. The numerical results 
of this example with sinc wavelet for k = 5 and N = 10, 20 given in [16] 
are reported in Table 1, which is comparable with our method. The value of 
the objective function for N = 2 is more accurate than its value for N = 3, 
because for N = 2, the quadratic functions of the state are approximated 
with better accuracy. The plot of approximate function u*(t) for N = 3 is 
given in Figure 1. The plots of the absolute error of state functions 7; and 
x2 for N = 3 are demonstrated in Figures 2 and 3. 


Table 1: Optimal value of J* for Example 2 


N 2 3 

J* | 8.55969 x 10-38 | 1.07179 x 10-*! 
N [16] 10 20 

if 6.88 x 10-2 6.88 x 10-2 


Figure 1: The approximate function of u* in Example 2 
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1 1 1 T 1 L 1 1 1 


-1.0 -0.5 r 0.5 1.0 


Figure 2: The absolute error of «} in Example 2 


-1.0 -0.5 0.5 1.0 


Figure 3: The absolute error of «3 in Example 2 


Example 3. Consider the following OCP governed by FIE systems 


subject to 
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The approximate and exact optimal control and state functions for N = 2 
are shown in Figure 4, 5 and 6, respectively. The optimal values of J* for 
different values of N are given in the Table 2. 


a 


1.0 

0.9 

0.8 

0.7 Exact function u(t) 

ont 2 (ikw.ULU™U™UUCU ti(<i‘OSOSCSCSC*C SS TS Approximate function u(t) 
0.5 

0.4 

4 4 4 4. r 4. 4 4. 4. 4 4 " 4. 4 4 4. 4. 4. 
0.0 0.2 0.4 0.6 0.8 1.0 


Exact function x;(t) 
----- Approximate function x; (t) 


0.0 0.2 0.4 0.6 08 1.0 


Figure 5: The exact and approximate optimal state 21(¢) for N = 2 in Example 3 
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Exact function x2(t) 


0.4 : ----- Approximate function x2(t) 


h 
7 
4 


0.0- 
ie 4 4 4 4 4 4 4 4 4 4 n 4 n 4 4. 4 4 n 4 did 
0.0 ny n no 4n 


Figure 6: The exact and approximate optimal state x2(t) for N = 2 in Example 3 


Table 2: Optimal value of J* for Example 3 


N 2 3 4 
J* | 1.10430 x 10-7 | 5.40995 x 10- | 1.11407 x 10-P 


Example 4 (Hanging chain). In this example, we investigate a problem 
from classical mechanics, which appears in power lines. If a flexible chain 
is hung from two points, it forms a U shape that is called a catenary curve. 
The resulting shape of the chain depends on the mass distribution along the 
chain. Consider the following FIE from [14], 


L 
eg | G(s, t)o(t)dt, (29) 


Here, x(s) and ¢(s) are demonstrated the displacement of the chain and the 
mass density of the chain respectively. g is the gravitational constant and 


st) 2a <t 
G(s, t) = Fie t Za Ze (30) 


In (30), T> is the constant tension. In this example, we are looking for 
the mass density distribution as the control function (t) that leads to a 
prescribed shape «(t) as its corresponding state. Therefore, the OCP problem 
is given as the minimization of the cost functional 


_ . g(tL? — t?) ‘ 2 
H(e,u)= [ (2 - 0) Le) =a da (31) 
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subject to equation (29). The optimal control and state solutions are 


o* (a) = a and a*(t) = de) In this example, we consider the following 
parameters as L = 1, g = 9.8, TO = 1. By applying the proposed algorithm 
for N = 3, the graphs of the exact and approximate controls and trajectories 
are given in figures 7 and 8. The optimal value of J* is 3.59596 x 1074. The 


graph of error for control and state function is given in figures 9a and 9b. 


—— Exact functon 


o.4f 4 oee-- Approximat2 functicn 


—— Exact function 


Figure 8: The exact and approximate optimal state in Example 4 
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(a) The error of state function in Example 4 
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(b) The errors of control function in Example 4 


Figure 9 


Example 5. Find the optimal control u* and corresponding optimal state 
x* that minimize objective function 


2 


io p= | (21 (€) —sin(t))? + (aro(t) — e-*)? + (u(t) — 2)? at 
0 


subject to 
xi(t) = h(t) — i (—t(u(s))?a1(s) + tu(s)29(s)) ds, 
X2(t) = h(t) — | (t(u(s) + 1)a1(s) + t?u(s)ao(s))ds. 
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. 23 5 2sin(5 107 
where hy(¢) = —sin(5t) 5 Ab + sm), 4 1Ory — ae 
e~ 3+ 4 2o08(5) ,_ sin(5),_£_ de" 424 "| The exact solutions are «7 (t) = sin(€), 


x(t) = e~* and u*(t) =t. Figures 10, 11 and 12 are the graph of exact and 
approximate solutions for N = 3 and N = 5, respectively. 


4 


eL 


1.0 


—_ 
0.0 


Exact function 


Figure 12: The accurate and proximate optimal state x2 in Example 5 


Iran. j. numer. anal. optim., Vol. 13, No. 2, 2023,pp 243-261 


Ebrahimzadeh 258 
Example 6. Find the optimal solutions that minimize objective function 
' 2 2 
min J = | (ai(t) — (t+ e"))” + (xa(t) — e’) 
0 
2 2 
+ (x3(t) — (1 + cos(t)))” + (u(t) — e')” dt 
subject to 
1 
x1(t) +f 3su(t)xo(s)ds =t + 4e’, 
0 
1 
£2(t) — i 6te°x1(s) — t?23(s)ds = e* — 3t(e? + 1) + t?(sin(1) +1) 
0 
lL: 
x3(t) — | (4(a3(s) — 1) + %2(s)ds = cos(t) — (e+ 4sin(1)) + 2 
0 
The exact solutions of this example are xj(t) = t+ e’, v3 (=e, 2) = 
1+ cos(t) and u*(t) =e’. The absolute errors of u*, xt, x} and x3 for N = 3, 
N =4and N = 5 are represented in Tables 3 and . 
Table 3: Absolute errors of control and state functions in Example 6 
N t=0 t=0.25 t=0.5 t=0.75 t=1 
3 | 7.51115 x 10-7 | 2.45764 x 10-7 | 4.67421 x 10-7 | 2.56842 x 10-7 | 9.41652 x 10-7 
Ex, | 4 | 3.99592 x 10-° | 7.55614 x 10~° | 5.88693 x 10~" | 9.47055 x 10° | 486114 x 10~° 
5 | 1.83038 x 10~© | 9.19475 x 1077 | 8.9541 x 10-7 | 9.96471 x 107-7 | 2.13043 x 107° 
3 | 4.91524 x 10-7 | 2.60616 x 10-7 | 5.32509 x 10-7 | 3.01636 x 10-7 | 5.72133 x 10-7 
Ex, | 4 | 2.49372 x 107° | 1.28883 x 107° | 1.78345 x 1077 | 1.4376 x 107° | 2.98331 x 107° 
5 | 1.05009 x 107© | 1.08341 x 107° | 1.13318 x 107% | 1.15422 x 107° | 1.19212 x 107° 
3 | 3.36834 x 10-7 | 1.04755 x 10-7 | 3.23772 x 10-7 | 9.70093 x 10~* | 3.05991 x 1074 
E,, | 4 | 6.42153 x 10~® | 3.47200 x 10~® | 9.92872 x 10~© | 4.26175 x 10~® | 8.97689 x 107° 
5 | 6.20439 x 1077 | 6.10935 x 1077 | 5.82297 x 107-7 | 5.87852 x 10-7 | 5.74393 x 1077 
3 | 4.34033 x 10-7 | 2.24118 x 10-7 | 6.13891 x 10-7 | 2.32091 x 10-7 | 5.65502 x 10-7 
E, | 4 | 1.93484 x 107° | 1.40343 x 107° | 7.62504 x 10~" | 1.68747 x 10~° | 2.46499 x 10-° 
5 | 8.13269 x 10-7 | 1.17233 x 10~® | 1.16398 x 10-® | 1.27319 x 10~® | 9.68081 x 107-7 


Table 4: Absolute errors of objective functions in Example 6 


N 3 A 5 
J* | 1.01824 x 10-7 | 2.36092 x 10-77 | 4.34198 x 10-78 


5 Conclusion 


The FIEs have rich physical and Engineering backgrounds. These equations 
got superior fondness across plenty of disciplines and broadly utilized in dy- 
namical machinery with chaotic or quasi-chaotic treatment. In this article, a 
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robust and efficient numerical approach based on Lagrange polynomials and 
collocation approach is applied to obtain the solutions of OCP governed by 
FIE system. At last, the prescribed problem is converted to a NLP. The re- 
sulted NLP is solved by optimization algorithm. The preciseness of Lagrange 
collocation approach can be readily concluded from the ameliorated results of 
our presented approach. The propounded approach can be extended for the 
approximation of solutions in OCPs involving FIE with fractional [5, 6, 12, 27] 
and singular [15] kernels, however some modifications will be needed in the 
method and convergence analysis. 
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